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SUMMARY 

The  contents  of  this  report  serve  to  fulfill  the  requirements  of  the  Phase  I  SBIR 
contract  #F49620-93-C-0021  between  AFOSR/PKC  and  C  &  C  Technologies,  Inc. 
entitled  "Integrated  Intelligent  Structures  for  Suppressing  Static  Aero  thermoelastic 
Deformations  and  Flutter  of  Panels".  The  ultimate  goal  of  this  project  is  to  develop  a 
commercial  package  of  methodology,  information,  and  criteria,  along  with  documented 
computer  programs,  for  designing  integrated  intelligent  panel  structures  to  meet  static 
aerothermoelastic  deformation  and  flutter  requirements.  To  accomplish  this  we  have,  in 
Phase  I,  investigated  the  properties  and  behaviors  of  integrated  intelligent  structures 
consisting  of  advanced  composite  laminates,  embedded  shape  memory  alloys  (SMA)  and 
piezoelectric  materials,  and  two  dynamic  control  systems.  In  this  Final  Technical  Report, 
the  results  of  the  Phase  I  effort  are  presented  in  detail  to  demonstrate  achievement  of  the 
objectives  of  the  project  and  to  show  how  this  innovative  integrated  intelligent  structure 
can  be  fully  exploited  in  Phase  II  and  beyond. 

By  successfully  meeting  or  exceeding  all  of  the  Phase  I  objectives,  we  have:  (1) 
Designed  optimally  shaped  piezoelectric  actuator  configurations  and  two  linear 
controllers  for  active  flutter  suppression;  (2)  Synthesized  a  simple  model-independent 
signature-based  controller  using  collocated  sensor/actuator  for  suppressing  nonlinear 
flutter  motion;  (3)  Developed  constitutive  relations  for  thermally  stressed  composite 
laminates  with  embedded  shape  memory  alloys  for  suppressing  static  aerothermoelastic 
deformations  and  passive  flutter  control;  (4)  Demonstrated  panel  flutter  suppression 
through  simulations  using  integrated  piezoelectric  active  control  and  SMA  passive 
control  for  laminated  composite  panels;  and  (5)  Included  large-deflection  nonlinearity  in 
the  control  simulations. 

In  Phase  I,  nonlinear  panel  flutter  suppression  has  been  studied  by  using  both  the 
finite  element  method  and  the  classical  continuum  method.  The  results  show  that  (1)  The 
bending  moment  is  more  effective  than  inplane  force  in  active  flutter  suppression  using 
piezoelectric  actuation.  (2)  With  a  pair  of  small  piezoelectric  layers  patched  at  the  leading 
edge,  the  critical  dynamic  pressure  can  be  increased  about  three  times  by  using  linear 
optimal  control  and  increased  about  two  times  by  using  the  model-independent  signature - 
based  control.  (3)  With  embedded  SMA  fibers,  the  critical  buckling  temperature  of  the 
composite  panel  is  increased  10  to  20  times  that  of  the  panel  with  no  SMA,  and  the  static 
aerothermoelastic  deformation  is  eliminated.  (4)  Flutter  stability  region  is  enlarged 
significantly  with  SMA.  (5)  Panel  flutter  suppression  for  integrated  composite  structures 
can  be  achieved  effectively  by  using  piezoelectric  materials  at  low  temperatures  and 
SMA  at  high  temperatures. 
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=  mcxlal  amplitude 
=  plate  bending  stiffness 
=  Young's  modulus 
=  Shear  modulus 
=  Mach  number 
=  number  of  the  actuator  sets 
=  moment  and  inplane  force  in  k-th  direction 
=  free  stream  airflow  speed 
=  positive  semi-definite  state  penalty  matrix 
=  positive  definite  control  penalty  matrix 
=  operation  temperature 
=  temperature  change 
=  Airy  stress  function 
=  nondimensional  modal  amplitude 
=  panel  length 
=  panel  width 

=  nondimensional  aerodynamic  damping 
=  mechanical  and  electric  coupling  coefficient 
=  electric  field 
=  panel  thickness 
=  aerodynamic  pressure 
=  time  and  nondimensional  time 
=  transverse  panel  deflection 
=  streamwise  coordinate 
=  spanwise  coordinate 
=  normal  coordinate 
=  Poisson’s  ratio 
=  density 
=  stresses 
=  strains 

=  nondimensional  dynamic  pressure 

=  system  aerodynamic  influence  matrices 
=  system  stiffness  matrices 
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=  system  mass  matrices 
=  system  aerodynamic  damping  matrices 
=  nondimensional  aerodynamic  damping  factor 
=  scale  factors  for  the  position  and  velocity  sensors 
=  position  and  velocity  feedback  gains 
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austenite  start  temperature 
austenite  finish  temperature 

=  the  first  and  second  nonlinear  stiffness  matrices 
panel  deflection 


Subscripts 


a  =  air 

b  =  bending 

cr  =  critical 

m  -  membrane 

max  =  maximum 

p  =  piezoelectric  material 

s  -  base  structure 

AT  =  temperature  change 


I.  INTRODUCTION 

Panel  flutter  at  supersonic  speeds  has  been  encountered  in  the  operation  of  aircraft  and 
missiles.  Panel  flutter  differs  from  the  conventional  lifting  surface  or  wing  flutter  in  at  least 
three  important  respects:  (1)  it  is  entirely  a  supersonic/hypersonic  phenomenon;  (2)  the 
airflow  is  only  on  one  side  of  the  panel;  and  (3)  the  structural  nonlinearities  due  to  large 
transverse  deflections  of  panels  tend  strongly  to  exhibit  large- amplitude  limit-cycle 
oscillations. 

The  earliest  reported  structural  failures  that  could  be  attributed  to  panel  flutter  were  the 
failure  of  the  60  to  70  early  German  V-2  rockets  during  World  War  II  [1,  2].  A  most 
recent  panel  flutter  failure  was  reported  by  Rick  Baker,  F-117A  Deputy  Chief  Engineer, 
Lockheed-Burbank  in  his  keynote  address  at  the  AIAA  Dynamics  Specialist  Conference 
[3].  After  the  flight  tests  of  the  F-l  17 A  Stealth  Fighter,  cracks  due  to  flutter  were  found  in 
almost  half  of  the  laminated  composite  skin  panels.  Those  panels  were  redesigned  and 
stiffened  and  thus  a  tremendous  weight  penalty  was  paid.  The  additional  weight  resulted  in 
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smaller  payloads  and  reduced  flight  speeds,  and  hindered  aircraft  maneuverability,  which  is 
extremely  important  for  fighter  aircraft  during  combat 

At  present  a  variety  of  high  speed  flight  vehicles  are  either  under  development  or  being 
considered  for  development  such  as  the  YF-22  Advanced  Tactical  Fighter  (ATF),  the 
National  Aero-Space  Plane  (NASP),  and  the  High  Speed  Civil  Transport  (HSCT).  All 
these  vehicles  feature  increased  performance  in  comparison  with  currently  existing  aircraft. 
Improvement  in  performance  characteristics  such  as  higher  speeds  (Mach  numbers),  greater 
maneuverability,  larger  payload  and  increased  efficiency  can  be  expected  to  lead  to  new 
panel  flutter  problems  arising  from  higher  dynamic  pressure  loads  (or  airflow  speed), 
higher  temperatures  and  large  deflections  of  light  weight  composite  structures.  When  a 
flight  vehicle  travels  at  high  speeds,  the  onset  of  flutter  may  be  imminent  but  usually  the 
dynamic  pressure  is  not  the  only  form  of  excitation.  An  aircraft  panel  is  subjected  to 
combined  aerodynamic  pressure  and  thermal  loads.  The  friction  from  the  surrounding  air 
causes  increased  heating  of  the  structure.  The  presence  of  high  temperature  loads  results  in 
flutter  motions  at  dynamic  pressures  much  lower  than  when  there  are  no  thermal  effects 
[4-8].  This  would  result  in  lower  flight  speed  limits.  In  addition,  the  temperature  rise  may 
also  cause  large  aerodynamic-thermal  deflections  of  the  skin  panels.  This  distortion  in 
shape  from  the  original  configuration  alters  the  lift-to-drag  characteristics  of  the  aircraft  and 
may  thus  lead  to  poor  performance. 

II.  PHASE  I  TECHNICAL  OBJECTIVES 

The  ultimate  goal  of  this  project  is  to  develop  a  commercial  package  of  methodology, 
information,  and  criteria,  along  with  documented  computer  programs,  for  designing 
integrated  intelligent  panel  structures  to  meet  static  aero  thermoelastic  deformation  and  flutter 
requirements.  To  accomplish  this  we  have,  in  Phase  I,  investigated  the  properties  and 
behaviors  of  integrated  intelligent  structures  consisting  of  advanced  composite  laminates, 
embedded  shape  memory  alloys  (SMA)  and  piezoelectric  materials,  and  two  dynamic 
control  systems.  Specific  objectives  include:  (1)  Initiate  panel  flutter  analysis  of  intelligent 
structures,  including  intelligent  actuator  design  using  piezoelectric  materials;  (2)  Study  the 
implementation  of  linear  optimal  controller  design  and  model-independent  signature-based 
controller  design;  (3)  Study  the  design  of  intelligent  actuators  using  shape  memory  alloys; 
(4)  Perform  integrated  design  and  development  of  actively-controlled  intelligent-structure 
panels  for  flutter  suppression. 


III.  ACTIVE  CONTROL  WITH  PIEZOELECTRIC  MATERIALS 
1.  Mathematical  Formulation 
1A.  Classical  Continuum  Method 

Piezoelectric  materials  can  develop  an  electrical  charge  when  subjected  to  a  mechanical 
strain.  The  converse  piezoelectric  effect,  the  development  of  mechanical  strain  when 
subjected  to  an  electrical  field,  can  be  utilized  to  actuate  a  structure  and  induce  inplane 
forces  and  bending  moments.  Thus,  actuation  of  a  structure  may  be  accomplished  at  the 
material  level.  The  basic  electric-mechanical  stain  and  stress  relation  is 


In  order  to  use  the  classical  method  to  investigate  the  panel  flutter  problem,  we  assume 
that  (a)  the  piezoelectric  material  used  here  is  isotropic;  (b)  the  piezoelectric  material  covers 
the  entire  top  and  bottom  surfaces  of  the  panel;  and  (c)  the  inplane  forces  induced  by  the 
piezoelectric  actuators  are  uniform.  The  piezoelectric  materials  will  produce  bending 
moments  and  inplane  forces  depending  upon  the  electric  field  applied  by  the  controller. 
These  inplane  forces  and  bending  moments  would  increase  the  critical  dynamic  pressure. 
Since  the  present  work  is  concerned  with  modification  of  the  structural  behavior  rather  than 
aerodynamic  precision,  the  first-order-piston-theory-type  relation  which  is  expressed  in  Eq. 
(2)  is  considered  to  provide  an  adequate  representation  of  the  supersonic  aerodynamic 
loads.  Consequently,  the  loads  are  obtained  from  the  relation 
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where  q  =  —  p„Vj  is  the  dynamic  pressure,  and  P  =  ^j M  2  - 1 .  The  nondimensional 

dynamic  pressure  is  given  by  A  =  — ^ — .  Another  commonly  used  nondimensional 

parameter  is  the  mass  ratio  which  is  defined  as  jU  =  where  mo  is  the  average  mass 

m0 

density  of  panel. 

The  piezoelectric  laminates  are  placed  symmetrically  with  respect  to  the  midplane,  so 


that  bending  and  stretching  deformations  are  decoupled.  We  also  assume  that  Poisson's 
ratio  is  the  same  for  both  materials,  i.e.  v,  =  \p  =  v  ,  consequently  the  nonlinear  equations 

of  motion  can  be  greatly  simplified  as, 
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where  the  definitions  of  variables  (see  Nomenclature)  follow  those  of  reference  [9].  For 
the  panel  flutter  analysis,  only  the  first  spanwise  mode  is  considered.  The  flutter  deflection 
of  a  simply  supported  plate  can  then  be  approximated  by 


?A-s,n(vM?)' 


(5) 


By  using  Galerkin’s  method,  the  governing  ordinary  nonlinear  differential  equation  in 
nondimensional  form  becomes  [9] 


+  C-j-  +  (K,  +  K  pUm  +  M  +  K„)Z  =  GU*. 
dir  dx 


(6) 


where  Z  is  the  vector  of  generalized  coordinates  an{=  /  h),  t  is  a  nondimensional  time, 
K qq  is  the  nonlinear  stiffness  matrix,  Um  and  Uj  are  normalized  control  variables  for 
inplane  force  and  bending  moment  respectively,  and  other  terms  are  straight  forward.  It  is 
noted  that 


c3m»x 


(7) 


where  e3m„  is  the  maximum  operating  electric  field  of  the  piezoelectric  material.  Because 
of  the  eZmilL  limitation,  for  a  given  Um,  there  is  a  constraint  for  the  normalized  control 

variable 
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IB.  Finite  Element  Method 


Using  the  finite  element  approach,  the  panel  can  be  studied  for  more  complicated 


geometries,  support  conditions  and  composite-laminate  arrangements.  The  nonlinear  finite 
element  equations  of  motion  are  derived  from  the  principle  of  virtual  work.  For  a  laminated 
aircraft  panel  consisting  of  fiber-reinforced  composite  and  piezoelectric  layers  subjected  to  a 
temperature  change  of  AT(x,y,z),  the  system  equations  of  motion  are  given  by 

-Um]{W}  +  —  [C]{W}  +  (A[A  ]  +  [K]  - [KNAT] - [KNe]  +  [Nl]  +  [N2]){W]  =  {PAX}  +  {Pe} 

K  0) o 

(9) 

where  the  definitions  of  symbols  (see  Nomenclature)  follow  those  of  reference  [10]. 
Because  of  the  large  number  of  degrees-of -freedom  usually  involved  with  finite  element 
model,  Eq.  (9)  is  not  suitable  for  control  to  suppress  the  panel  deflection  {W}.  A  modal- 
coordinate  transformation,  therefore,  is  used  to  convert  it  to  a  reduced  number  of  modal 
coordinates.  The  resulting  reduced  second-order  equation  is  shown  in  Eq.  (6).  Numerical 
results  for  nonlinear  panel  flutter  motions  obtained  from  the  classical  continuum  method 
previously  described  validate  results  obtained  from  the  finite  element  method. 

2.  Piezoelectric  Actuator  Configuration  Design 

A  typical  laminated  aluminum  panel,  shown  in  Figure  1,  contains  two  piezoelectric 
layers  of  equal  thickness  placed  symmetrically  with  respect  to  the  midplane  of  the  panel. 
The  actuator  set  can  be  stimulated  such  that  one  laminate  contracts  and  the  other  expands  to 
create  bending  moments  in  the  structure,  such  that  both  layers  exhibit  the  same  strain 
resulting  in  nonzero  inplane  forces,  or  a  combination  of  the  two.  Both  induced  bending 
moments  and  inplane  tensile  forces  can  increase  the  stiffness  of  the  controlled  structure. 


Figure  2  shows  three  piezoelectric  actuator  arrangements  considered  in  Phase  I.  They  are: 
one-set  patched,  two-set  patched,  and  one-set  small  patched  embedded  piezoelectric 
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laminates.  In  the  first  two  configurations,  an  isotropic  panel  is  completely  covered  by  the 
piezoelectric  layers.  For  the  two-set  piezoelectric  actuator  configuration,  each  set  actuator  can 
be  activated  by  different  controller.  Note  that  the  one-set  configuration  is  a  limiting  case  of 
the  two-set  configuration.  In  the  last  configuration  only  a  small  part  of  the  top  and  bottom 
layers  of  an  eight-layer  composite  panel  is  embedded  with  piezoelectric  materials  adjacent  to 
the  leading  edge. 


1-set  Piezo 
Isotropic  Panel 


2-set  Piezo 
Isotropic  Panel 


1-set  Small  Piezo 
Composite  Panel 


Panel  1  Panel  2  Panel  3 


Fig.  2  Piezoelectric  Actuator  Configuration  Design. 


3.  Linear  Control  Designs 
3 A.  Linear  Optimal  Control 

For  the  suppression  of  panel  flutter,  an  optimal  control  approach  based  on  the  linear 
optimal  control  theory  is  considered  first.  The  control  design  is  based  on  the  linear  model 
obtained  by  ignoring  the  nonlinear  term  of  Eq.  (6).  For  optimal  control  using  both  inplane 
forces  and  bending  moments,  the  stiffness  matrix  varies  with  U m  .  The  problem  can  be 
simplified  by  first  assuming  Um  to  be  a  constant.  The  control  variable  then  becomes  UA 
which  induces  bending  moment  only.  From  Eq.  (6),  the  first  order  dynamic  equation 
becomes 

—  =  AX  +  BU,  (10) 

dx 

where 

I"Z1  T  0  1  1  r  0  1 

X=  H  ’  A=  -M~'(K,  +  KpUm  +  KaX)  -M  'C  ’  B  LM-'GJ’  U-Ut’ 


The  linear  quadratic  performance  index  for  optimal  control  can  be  formulated  as 


J  =  -  J“(XrQX  +  UrRU)dr,  (11) 

where  Q  is  a  positive  semi-definite  state  penalty  matrix  and  R  is  a  positive  definite  control 
penalty  matrix.  From  the  optimal  control  theory,  the  optimal  controller  for  this  linear 
quadratic  regulator  problem  is  obtained  by 

U  =  -R_1BrPX,  (12) 

where  P  is  a  positive  definite  matrix  obtained  from  the  following  Riccati  equation  [10] 

AtP  +  PA-PBR~1BtP  =  -Q.  (13) 


The  weighting  matrix  R  is  chosen  as  an  identity  matrix  times  a  positive  constant,  and  the 
weighting  matrix  Q  is  chosen  by  the  energy  weighting  method  in  this  study 


Q  = 


(14) 


For  a  different  constant  Un,  an  optimal  feedback  control  gain  for  Ub  can  be  obtained. 

Through  numerical  simulations  using  nonlinear  equations  of  the  panel,  the  critical  dynamic 
pressure  can  be  obtained.  Finally,  an  optimal  set  of  Um  and  Ub  is  obtained  by  comparing 
the  results  obtained  from  different  constant  Um. 

3B.  Model-Independent  Signature-Based  Control 


In  this  section,  we  study  the  feasibility  of  the  model-independent  signature-based 
control  design.  For  linear  optimal  control,  we  need  full  state  feedback  and  an  accurate 
model  to  design  the  controller.  In  general,  it  is  very  difficult  to  design  a  set  of  sensors  to 
measure  all  the  states  which  include  the  position  and  velocity  of  each  mode.  If  the  full  state 
information  can  not  be  directly  obtained  from  sensors,  we  need  a  state  estimator  to  estimate 
all  the  states  indirectly.  On  the  other  hand,  an  accurate  model  is  also  very  difficult  to 
obtain.  Either  the  state  estimation  error  or  the  modeling  error  may  degrade  the  performance 
of  the  controller.  This  motivates  us  to  design  the  model-independent  signature-based 
controller  to  avoid  these  two  problems.  The  term  model-independent  means  that  the 
feedback  gains  in  the  controller  design  can  be  determined  without  any  model  information. 
Output  feedback  is  used  to  replace  the  state  feedback.  The  term  signature-based  means  that 
the  model  is  still  important  to  decide  the  optimal  locations  for  collocated  sensors/actuators. 


This  controller  design  is  very  similar  to  the  passive  damping  and  stiffness  design  for 
structural  systems  [11].  As  we  know,  when  a  mass-spring-dashpot  is  attached  to  any 
mechanical  system,  the  damping  and  stiffness  of  the  system  are  almost  always  augmented 
and  the  stability  of  the  system  is  thus  improved.  The  question  arises  as  to  whether  there  are 
any  feedback  control  designs  using  sensors  and  actuators  which  behave  like  the  passive 
mass-spring-dashpot.  These  kinds  of  controllers  sometimes  are  called  virtual  passive 
controllers  and  can  be  achieved  by  using  collocated  sensors/actuators  and  direct  negative 
feedback  controls.  Recall  a  simplified  linear  model  equation  from  Eq.  (6) 

MZ  +  CZ  +  KZ  =  GU.  (15) 

For  collocated  sensors/actuators  and  a  direct  negative  feedback  control,  one  has 
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U  =  -[*/>  ^,]Y  =  -(/AGrZ  +  /AGrZ).  (16) 

where  f  p  and  /v  are  scale  factors  for  the  position  and  velocity  sensors,  respectively,  and 
kp  and  kv  are  position  and  velocity  feedback  gains,  respectively.  The  closed-loop  linear 

system  becomes 

MZ  +  (C  +  /AGt)Z  +  (K  +  fpkp  GT)Z  =  0.  (17) 

The  augmented  damping  fykvGT  and  stiffness  fpkp Gr  depend  on  fcv,  kp  and  G.  Larger 
kv  and  kp  can  provide  larger  damping  and  stiffness  but  need  larger  control  power. 

Meanwhile,  G  depends  on  the  size,  shape  and  location  of  the  collocated  sensors/  actuators. 
To  achieve  proper  augmented  damping  and  stiffness  for  suppressing  the  panel  flutter,  we 
have  to  find  the  best  G  by  investigating  all  possible  configurations  of  the  collocated 
sensors/actuators.  The  knowledge  of  the  model  can  be  used  for  the  optimal 
sensors/actuators  configuration  design. 

4.  Numerical  Results 

4A.  Results  for  Isotropic  Panel 

In  this  section,  numerical  results  are  presented  for  suppressing  flutter  of  an  isotropic 
panel  using  the  linear  optimal  control  and  the  model-independent  signature-based  control. 
The  panel  studied  is  a  simply  supported  square  aluminum  panel  with  two  piezoelectric 


laminates  which  cover  the  two  surfaces  of  the  panel  (Panel  1  of  Fig.  2).  The  physical 
parameters  and  the  geometry  of  the  panel  are  shown  in  Appendix. 

First,  optimal  control  design  is  performed.  Figure  3  shows  time  response  of  the 
maximum  panel  deflection  and  the  feedback  control  effort.  The  nondimensional  dynamic 
pressure  X  is  1400,  and  the  inplane  force  is  set  to  zero.  Time  response  of  the  panel 
deflection  shows  a  limit  cycle  motion  at  the  beginning.  After  the  controller  is  activated,  the 
panel  deflection  is  quickly  suppressed. 

The  inplane  force  can  change  the  panel  stiffness.  Inplane  tension  increases  the 
stiffness,  and  inplane  compression  decreases  the  stiffness.  Figure  4  shows  the  inplane 
force  effect  on  the  Xcr  with  or  without  bending  moment  used  for  optimal  feedback  control. 

The  range  of  control  authority  for  bending  moment  varies  with  the  inplane  force.  When  the 
inplane  force  is  set  to  the  maximum  or  minimum  value,  no  bending  moment  can  be  induced 
because  of  the  limitation  shown  in  Eq.  (8).  As  the  inplane  force  approaches  zero,  the  Xcr 

increases  dramatically,  and  the  flutter  free  region  is  enlarged.  With  zero  control  of  inplane 
stress  and  maximum  control  of  bending  moment,  the  Xcr  can  be  increased  to  three  times  the 

value  for  the  uncontrolled  panel.  Control  of  the  bending  moment  is  obviously  the  more 
effective  way  to  suppress  flutter. 

The  one-set  patched  actuator  is  an  extreme  case  of  the  two-set  patched  actuator  (Panel  2 
in  Fig.  2).  If  normalized  separating  position  is  set  to  0  or  1,  then  the  two-set  patched 
actuator  becomes  a  one-set  patched  actuator.  Figure  5  shows  the  variation  of  Xcr  with  the 

normalized  separation  position  of  a  two-set  patched  actuator  and  the  control  weighting 
factor  r.  The  inplane  force  is  set  to  zero.  The  result  shows  that  the  two-set  actuator  can 
perform  better  than  a  one-set  actuator.  Table  1  also  shows  that  if  only  the  first  set  is 
activated  (x/a  =  0.3,  for  example),  the  still  can  be  increased  more  than  three  times  the 

value  514  for  the  uncontrolled  panel.  In  fact,  one-set  activated  with  x/a  as  low  as  0.3  is 
more  effective  than  with  x/a  =  1.0  (i.e.,  equivalent  to  Panel  1). 

Optimal  location  for  small  patched  actuators  is  also  studied  by  analyzing  the  norm  of  the 
optimal  feedback  gain  at  various  locations  of  the  panel  (see  Fig.  6).  The  finite  element 
model  is  meshed  by  10x6  rectangular  elements.  The  nondimensional  dynamic  pressure  is 
1000.  Numerical  results  show  that  if  the  actuator  is  patched  at  the  locations  with  higher 
norm,  it  can  suppress  flutter  more  effectively.  If  the  piezoelectric  materials  are  patched  at 
the  locations  with  the  norms  higher  than  120  (adjacent  to  the  leading  edge),  the  critical 
dynamic  pressure  can  be  increased  about  four  times  of  the  value  without  piezoelectric 
actuation. 
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Fig.  3  Time  History  of  Panel  Deflection  and  Control  Effort 


Fig.  4  Critical  Dynamic  Pressure  Xcr  with/without  Bending  Moment  Control. 
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Fig.  5  Critical  Dynamic  Pressure  vs.  Normalized  Separating  Position  (x/a)  and 
Control  Weighting  Factor  r  for  a  Two-Set  Patched  Actuator. 


Table  1  Comparison  of  Xcr  for  2-Set  Piezo,  Isotropic  Panel,  with  R  =  10,000  x  I 
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Fig.  6.  Norm  of  the  Optimal  Feedback  Gain  for  Upper  Half  Panel 


Fig.  7.  Optimal  Collocated  Sensor/Actuator  Configuration  (Shaded  Part)  for  Model- 
Independent  Signature-Based  Control. 

Next,  the  model-independent  signature-based  control  design  is  studied.  The  finite 
element  model  is  meshed  by  10x6  rectangular  elements.  Optimal  sensor/actuator  locations 
(see  Fig.  7)  are  obtained  by  investigating  all  possible  configurations  of  the  collocated 
sensors/actuators.  The  feedback  gains  are  not  only  model  independent  but  also 
aerodynamic  independent.  The  result  shows  that  the  nondimensional  critical  dynamic 
pressure  can  be  increased  from  514  up  to  about  1 100. 

Since  this  model-independent  controller  does  not  require  the  knowledge  of  the  model 
and  the  state  information,  it  is  very  simple  and  attractive  for  experimental  validation  in 
Phase  II.  However,  the  design  of  the  size,  shape  and  location  of  the  collocated  sensors/ 
actuators  to  maximize  the  critical  dynamic  pressure  needs  further  study  by  using  the 
knowledge  of  the  model.  A  new  finite  element  analysis  based  on  triangular  elements 
instead  of  rectangular  also  needs  to  be  developed  to  further  improve  the  sensor/actuator 


configuration. 


4B.  Results  for  Composite  Panel 

A  simply  supported  square  thin  composite  panel  is  studied.  The  panel  is  of  dimension 
12"xl2"x0.05",  and  is  made  up  of  eight  layers  with  lamination  angle  [0/45/-45/90],.  All 
inplane  displacements  are  considered  to  be  zero  at  the  edges  (u=v=0).  The  material 
properties  for  the  panel  are  shown  in  Appendix.  A  one-set  small  patched  piezoelectric  layer 
is  embedded  within  the  composite  top  and  bottom  0-deg.  layers  (see  Fig.  2).  These 
piezoelectric  layers  have  an  isotropic  behavior.  Their  properties  are  also  shown  in 
Appendix. 

The  panel  is  modeled  by  a  10x6  array  of  C1  conforming  rectangular  bending  element. 
The  modal  equation  of  motion,  Eq.  (6),  is  solved  by  Runge-Kutta  method.  The  inplane 
forces  induced  by  the  piezoelectric  layers  are  null,  since  the  voltages  applied  on  top  and 
bottom  piezoelectric  patches  are  with  the  same  values  but  different  signs.  Figure  8  plots  the 
maximum  panel  deflection  vs.  nondimensional  dynamic  pressure  for  different  numbers  of 
modes  retained  in  Eq.  (6).  The  six-mode  representation  appears  to  give  a  converged  limit- 
cycle  result. 


Fig.  8  Maximum  Panel  Deflection  with  Fig.  9  Maximum  Panel  Deflection  with 
Different  Number  of  Modes.  Different  Piezo  Sizes. 


The  effect  of  the  embedded  piezoelectric  actuator  size  has  also  been  studied  (Panel  3  in 
Fig.  2),  and  the  results  are  illustrated  in  Fig.  9.  The  larger  the  inactive  piezo  patch,  the 
lower  the  dynamic  pressure  is,  because  inclusion  of  piezoelectric  materials  reduces  the 
panel  stiffness.  However,  if  piezoelectric  actuator  is  activated  using  linear  optimal  control, 


the  critical  dynamic  pressure  can  be  increased  from  279  to  686  for  the  2.4"x8"  piezo  size 
and  from  298  to  429  for  the  1.2"x4"  piezo  size.  The  optimal  shape  and  size  of  the 
piezoelectric  actuator  will  be  studied  in  Phase  II. 


IV.  PASSIVE  CONTROL  WITH  SMA 

It  is  known  that  the  compressive  thermal  stresses  induced  by  high  temperature  generally 
reduce  panel  stiffness  and  hence  reduce  the  critical  dynamic  pressure.  Lower  speed-limit  of 
the  aircraft  could  result.  This  is  illustrated  in  Fig.  10  for  the  same  laminated  panel 
discussed  in  the  previous  section  but  without  piezoelectric  material.  Shape  memory  alloys 
(SMA),  however,  have  three  interesting  features  in  a  thermal  environment.  When  an 
initially  strained  elongated  SMA  fiber  is  heated  above  its  characteristic  transformation 
(austenite  finish  Af)  temperature,  it  will:  (1)  return  to  its  original  length  or  produce  a 
recovery  tensile  stress  in  the  fiber  with  an  edge-immovable  (strain  restrained)  boundary;  (2) 
increase  its  Young's  modulus  by  a  factor  of  three  or  four;  and  (3)  greatly  increase  its 
yielding  strength.  These  characteristics  are  shown  in  Fig.  11  [12]  for  55-Nitinol  alloy. 

The  aerodynamic  heating  exerted  on  the  skin  panels  of  a  supersonic  vehicle  will  cause 
panel  temperatures  well  above  the  austenite  finish  temperature  Af.  Thus  the  three  features 

of  the  SMA  will  be  ideal  for  passive  control  of  flutter  by  embedding  the  SMA  fibers  in  a 
laminated  composite  panel.  The  SMA  will  reduce  the  adversary  thermal  stress  effects  on 
panel  flutter  and  increase  the  critical  dynamic  pressure  of  the  panel. 

1.  Mathematical  Formulation 

In  previous  studies  of  SMA  reinforced  composites  [14,15],  the  thermal  stress  effect  in 
the  composite  matrix  was  neglected.  For  the  application  of  SMA  in  a  thermal  environment, 
this  effect  has  to  be  considered  in  the  formulation.  The  constitutive  equation  for  a 
composite  laminate  with  embedded  SMA  fiber,  which  includes  the  thermal  effects,  is 
developed  in  this  study. 

1A.  Properties  of  a  Composite  Lamina  with  Embedded  SMA  Fibers 

The  stress-strain  relations  of  a  composite  lamina  (for  example,  graphite-epoxy)  with 
embedded  SMA  fibers  in  material  principal  directions  are  [16] 
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Fig.  10  Limit-Cycle  Amplitude  of  a  Thermally  Loaded,  Simply  Supported,  Laminated 
Panel. 


Fig.  11  Shape  Memory  Alloy  Material  Property. 
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where  the  reduced  stiffness  matrix  [<2m]  is  evaluated  by  EXm,  E2m,  G12m,  and  v12m  of  the 
composite  matrix,  and  Vn  and  Vs  are  volume  fractions  of  composite  matrix  and  SMA 
fibers.  The  reduced  stiffness  matrix  [£)]’  is  evaluated  by  Ev  E2,  G12,and  v12oftheSMA 

reinforced  composite  [16], 


a2=Vma2m, 


( p  r  \_  {^2m^s^\2mPs) 

[  [(£2».g,2.)v,+(£;.g;)v.]' 


(20) 


The  superscript  (*)  indicates  that  the  quantity  is  temperature  dependent,  and  the  recovery 
stress  <r*  of  SMA  depends  on  the  temperature  as  well  as  the  initial  strain  er  (Fig.  1 1). 

For  a  general  k-th  layer  of  lamination  angle  0,  the  stress-strain  relations  in  laminate 
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The  resultant  stress  are  defined  as 

{N,M}  =  J_“{a}k(U)dz. 


(23) 
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IB.  Finite  Element  Equations  of  Motion 


By  using  the  principle  of  virtual  work,  the  governing  equations  of  motion  are  derived 
for  the  panel  flutter  with  the  effects  of  thermal  and  SMA  recovery  stresses  as, 


-4[M](tV) + %G](H')  +  aiA.] + m 
<i>0 

+[Af„] + |[W1]  +  |[iV2]){W'}  =  {Par}  -  {?„},  (24) 


where  [G]  and  [AJ  are  the  aerodynamic  damping  and  influence  matrices,  [KNAT]  and  [Ka] 

are  the  thermal  geometric  and  SMA  recovery  stress  stiffness  matrices,  [Nl]  and  [N2]  are 
the  first  and  second  order  nonlinear  stiffness  matrices,  and  {PAr}  and  are  the  thermal 

load  and  SMA  recovery  force  vectors.  The  nondimensional  dynamic  pressure  "K  as  defined 
from  Eq.  (2), 


A  = 


2  qa3 


(25) 


where  D  =  Dno  is  the  Dn  of  the  composite  material  (Vm=l)  with  all  fibers  aligned  in  0°  or 
x-direction.  The  aerodynamic  damping  from  Eq.  (2)  is 


8a  ph(00{Ml- 1)3/2 


(26) 


1C.  Solutions 


The  equation  of  motion,  Eq.  (24),  is  a  set  of  ordinary  differential  equations  with  respect 
to  time  t.  However,  the  force  vectors  {^Ar}  and  are  time  independent.  Thus  the 

total  solution  is 

{W}  =  {W}s+{W}t,  (27) 

where  {W}s  is  the  time-independent  particular  solution  which  yields  the  aerothermal 
postbuckling  deflection  with  SMA  recovery  stress  effect,  whereas  {W}t  is  the  time- 
dependent  homogeneous  solution  of  linear  flutter  motions.  Substituting  Eq.  (27)  into  Eq. 
(24),  the  following  two  equations  are  obtained: 

(A[AJ  +  [ K ]  -  [Knat]  +  [Ka]  +  i[M],  +  i[N2],  ){W}S  =  {PAT}  -  {Pa}  (28) 

-\[M]{W},  +  ^-[G]{W}t  +  (X[Aa)  +  [K] - [Knat]  +  [Ka]  +  [Nl],  +  [N2]s)[W}t  =  0 .  (29) 

°>o  <»o 


The  Newton-Raphson  iteration  method  is  well  suited  to  determine  {W)s  from  Eq.  (28), 
and  the  coalescence  of  the  lowest  two  modes  determines  the  critical  dynamic  pressure  from 
Eq.  (29).  Detailed  solution  procedures  are  referred  to  Refs.  [7,  8]. 

2.  Numerical  Results  and  Discussions 

A  simply  supported  SMA  fiber  embedded  [0/45/-45/90]  s  graphite-epoxy  square  panel 
(a=b=12  in,  h=0.05  in)  is  studied.  The  SMA  fibers  are  homogeneous  throughout  the  panel 
in  the  same  directions  as  of  the  graphite  fibers.  The  material  properties  of  the  graphite- 

epoxy  and  SMA  fiber  are  shown  in  Appendix. 

It  is  known  that  the  recovery  stress  a*  (see  Fig.  11)  is  a  nonlinear  function  of  the 

temperature  and  the  initial  strain,  at  high  temperatures  it  remains  almost  unchanged.  At  180 
°F,  the  recovery  stress  can  be  calculated  from  the  expression  [17] 

<r*=  5.34xl06er3  -  4.89xl05er2+  1.68x10^  (MPa),  er<5%.  (30) 

In  the  finite-element  results,  it  is  assumed  that  the  recovery  stress  is  constant  and  can  be 
calculated  from  Eq.  (30)  when  T>Af. 

2A.  Critical  Buckling  Temperature  and  Thermal  Postbuckling 

Tables  2  shows  the  critical  buckling  temperatures  and  the  maximum  thermal 
postbuckling  deflections  of  the  [0/45/-45/90]s  graphite-epoxy  composite  panel  embedded 
with  SMA  fibers  of  various  initial  strains  and  volume  fractions.  It  is  observed  that  for 
T>Af  the  buckling  temperatures  increase  with  increasing  volume  fraction  and  with 
increasing  initial  strain  of  the  SMA  fibers.  At  low  temperatures  (T<As)>  the  critical 
temperature  is  decreased  with  the  increasing  of  SMA  volume  fraction,  because  Es  <  Elm. 

Large  thermal  postbuckling  deflection  could  alter  the  aerodynamic  flow  field  and  thus 
reduce  the  operation  performance  and  efficiency  of  the  aircraft  as  well  as  significantly 
increase  surface  temperatures.  The  elimination  or  reduction  of  the  thermal  postbuckling 
deflection  by  embedded  SMA  fibers  will  maintain  the  original  optimal  aerodynamic 
performance  designed  for  the  aircraft 

Figures  12  and  13  show  that  the  flutter  stability  region  is  enlarged  significantly  at  high 
temperatures  (T>Af)  with  the  increase  of  SMA  fiber  volume  fraction  and  initial  strain. 
Figure  14  shows  that  the  stability  region  decreases  slightly  at  low  temperatures  (T<AS) 
when  the  panel  has  more  SMA  fibers  because  E*  <  Elm  as  previously  noted. 


Table  2  Critical  Buckling  Temperatures  ATCT  (°F) 


SMA  Vol. 
Vs  (%) 

T<  As 

T>  Af 

Er  =  2% 

£j  =  3% 

ii 

ctf 

0 

21.4 

- 

- 

- 

10 

18.7 

152 

169 

185 

20 

16.3 

314 

354 

390 

30 

14.3 

523 

591 

653 

In  the  transient  temperature  range  (As<T<Af),  only  part  of  the  total  initial  strain  can  be 
recovered,  and  the  recovery  stress  can  not  be  estimated  from  Fig.  11.  Detailed  calculations 
showed  that  the  critical  dynamic  pressure  increases  linearly  with  the  temperature,  and  the 
stability  boundary  is  a  straight  line  between  (As,  Xs)  and  (Af,  Xf).  Figure  15  shows  the 
complete  flutter  stability  boundaries  for  the  composite  panel  studied  with  SMA  (Vs=0.2, 
Er=0.03)  and  without  SMA  (Vs=0).  This  certainly  demonstrates  that  SMA  can  greatly 
enlarge  the  stable  region  especially  at  high  temperatures. 

V.  INTEGRATED  INTELLIGENT  STRUCTURE 

In  sections  HI  and  IV,  the  use  of  piezoelectric  material  and  shape  memory  alloy  (SMA) 
to  suppress  panel  flutter  motions  have  been  investigated,  respectively.  The  dynamic 
equations  of  motion  have  been  derived.  The  primary  results  for  isotropic  as  well  as 
composite  panels  are  obtained.  If  limit-cycle  flutter  motions  occur,  they  can  be  suppressed 
quickly  within  the  critical  dynamic  pressure  with  piezoelectric  actuators  (Fig.  3).  The 
flutter  boundary  is  increased  tremendously  at  high  temperatures  with  the  use  of  SMA. 
Therefore,  it  is  suggested  that  the  piezoelectric  materials  are  to  be  used  to  actively  suppress 
the  panel  flutter  at  the  low  temperatures  (T<AS)  since  the  SMA  will  not  be  activated  in  this 
region.  On  the  other  hand,  the  SMA  should  be  adopted  to  enlarge  the  stability  flight 
envelope  at  the  high  temperatures  (T>Af). 

In  this  section,  both  piezoelectric  material  and  SMA  are  embedded  in  the  same 
laminated  composite  panel.  SMA  fibers  are  aligned  in  the  same  direction  as  the  graphite 
fibers  in  each  composite  lamina.  The  finite  element  time  domain  approach  is  adopted  in  this 
study. 

Finally,  consider  the  example  of  a  panel  with  only  20%  SMA  fibers  embedded  in  the 
composite  laminate  and  residual  strain  taken  to  be  3%,  and  with  embedded  piezo  of 
dimension  2.4"x8"  at  the  leading  edge.  The  critical  dynamic  pressures  are  shown  in  Fig. 
16  for  the  panel  at  various  temperatures.  By  using  integrated  SMA  and  piezo,  the  stability 


Fig.  12  Flutter  stability  boundaries  of  SMA/composite  plates  (£^=0.02)  with  different  SMA 
volume  fractions 


Fig.  13  Flutter  stability  boundaries  of  SMA/composite  plates  (Vs=0.3)  with  different  initial 
recovery  strains 


Fig.  14  Flutter  stability  boundaries  of  SMA/composite  plates  with  different  SMA  volume 
fractions  at  low  temperatures  (T<AS) 


Fig.  15  Complete  stability  boundary  for  SMA/composite  panel  with  Er=0.03. 
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Fig.  16  Complete  stability  boundary  for  a  composite  panel  with  or  without  piezo  and 
SMA(Vs=0.2  and  £^0.03). 

boundary  is  enlarged  throughout  the  temperature  range.  The  critical  dynamic  pressure  is 
therefore  increased.  When  T<AS,  the  SMA  is  not  active.  On  the  other  hand,  when  T>AS, 
SMA  is  automatically  activated,  but  with  relatively  small  recovery  stress.  When  T>Af,  the 
large  recovery  stresses  will  stiffen  the  panel  so  that  the  critical  dynamic  pressure  is  further 
increased.  In  this  region,  the  passive  control  using  SMA  is  automatically  achieved.  In  this 
study,  however,  the  piezoelectric  materials  are  assumed  to  be  independent  of  temperatures. 
The  effect  of  temperature  on  piezoelectric  material  will  be  included  in  Phase  H 

VI.  CONCLUSION 

The  ultimate  goal  of  this  project  is  to  develop  a  commercial  package  of  methodology, 
information,  and  criteria,  along  with  documented  computer  programs,  for  designing 
integrated  intelligent  panel  structures  to  meet  static  aerothermoelastic  deformation  and  flutter 
requirements.  To  accomplish  this  we  have,  in  Phase  I,  investigated  the  properties  and 
behaviors  of  integrated  intelligent  structures  consisting  of  advanced  composite  laminates, 
embedded  shape  memory  alloys  (SMA)  and  piezoelectric  materials,  and  two  dynamic 
control  systems. 

The  fundamental  conclusion  of  this  Phase  I  effort  is  that  it  is  possible  to  use  existing 
technology  to  develop  the  proposed  innovative  integrated  intelligent  structure  to  meet  static 
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aerothermoelastic  deformation  and  flutter  requirements  and  to  systemize  the  process  of 
designing  such  structures  to  meet  specified  requirement.  Panel  flutter  suppression  can  be 
achieved  for  integrated  intelligent  composite  structure  by  using  piezoelectric  materials  at 
low  temperatures  and  SMA  at  high  temperatures.  Optimal  shaped  piezoelectric  actuator 
configuration  design  shows  that  actuating  a  small  piezoelectric  layer  patched  adjacent  to  the 
leading  edge  can  suppress  the  flutter  effectively.  The  model-independent  signature-based 
controller  using  collocated  sensor/actuator  is  very  simple  and  easy  to  implement  as  an 
effective  active  controller.  With  embedded  SMA  fibers,  the  critical  buckling  temperature  of 
the  composite  panel  is  increased  10  to  20  times  as  compared  to  panel  with  no  SMA  and  the 
flutter  stability  region  can  be  enlarged  significantly. 

Under  Phase  I,  we  have  demonstrated  the  feasibility  of  all  the  critical  design  and 
nonlinear  analysis  of  the  innovative  intelligent  structure  for  integration  in  Phase  II  into 
designing  of  a  realistic  aircraft  intelligent  panel  for  ultimate  laboratory  and  wind  tunnel 
testings.  The  extensive  experience  of  methodology,  design  and  analysis  we  have  available 
will  enable  us  to  develop  integrated  intelligent  panel  structures  to  meet  static 
aerothermoelastic  deformation  and  flutter  requirements  and  respond  to  future  requirements. 

By  successfully  meeting  and  exceeding  all  of  the  Phase  I  objectives,  we  have 
demonstrated  an  innovative  solution  to  an  existing  need  for  suppressing  static 
aerothermoelastic  deformations  and  flutter  of  panels.  In  this  way.  Phase  I  has  served  its 
purpose  as  a  design,  analysis  and  demonstration  period  and  paved  the  way  for  complete 
laboratory  and  wind  tunnel  testing  in  Phase  II.  The  far  reaching  implications  of  the 
intelligent  structure  introduce  new  possibilities  in  both  the  government  and  commercial 
arenas.  It  is  these  areas  which  will  allow  widespread  development  and  dissemination  of  the 
intelligent  structure  and  a  commercial  package  for  its  design  under  Phase  HI  and  hence 
fulfill  the  goals  of  the  SBIR  program. 
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APPENDIX 


Properties  of  Panel  Used  in  Calculations 


Isotropic  Panel  (see  Figure  11 

a  =  12  in 

b  =  12  in 

h  =  0.08  in 

hb  =  0.05  in 

Eb=  1.0x10^  psi 

Vb  =  0.3 

pb=  0.2588x10-3  lb-s2/in4 

GraphiteZEpQxy  Panel 
Elm  =  22.5E+6  psi 
E2m  =  1.17E+6  psi 
GUm  =  0.66E+6  psi 

vi2m  =  0-22 

pm=  0.056  lb/in3 
aln  =  -0.04E-6  /°F 
a2m  =  16.7E-6  /oF 


piezoelectric  Material 
ap=3.0xl0‘6/eF 
dj;=dj2  =-7.48xl0'9  in/v 
e3max=  1.52X104  v/in 
hp  -  0.015  in 
£’p=9.0xl06  psi 
Vp  =  0.3 

pp=  0.7101x10-3  lb-s2/in4 
SMA-Nitinol 

E]  =  3.996E+6  psi,  (T<As) 

E\  =  11.994E+6  psi,  (T>Af) 

G,*  =  3.604E+6  psi,  (T<As) 

G;  =3.712E+6psi,  (T>Af) 

ps  =  0.2325  lb/in3,  v,  =  0.3 
As  =  100  °F,  Af  =  145  °F 
as  =  5.7E-6  /°F,  er  =  2~4% 


